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Abstract 

In this paper, we use Martin and Coulomb-Linear potentials and solve Lippman- 
Schwinger equation and then identify be energy levels. Moreover, we predict results 
for such energy levels as that of tt (in its short half-life) which is not observed. We 
showed our results are consistent with previous findings in literature. Also inves¬ 
tigating spectrum of eigen-values, we obtain stability interval for Yukawa-Linear 
potential. 

keyword: 

Schwinger equation. Mass spectrum. Strong potential models. Cornel potential, heavy 
quark 


*nt_physics @ yahoo.com 
1 monem@kashanu. ac. ir 
1 hadizade@ift. unesp .hr 


1 



1 Introduction 


Hadrons consist of quarks and antiquarks and follow the principles of Quantum Chorody- 
namics (QCD). Perturbative measuring of physical properties of hadrons does not always 
yield consistent results except for the deep inelastic region in which the coupling constant 
is influenced by asymptotic freedom and is weakened. Therefore, such theories as lattice 
gauge theory have emerged to consider non-perturbative effects. Among hadrons, mesons 
are favorable means of exploring strong interaction in the non-perturbative method with 
strong coupling. Hence, studying meson spectrum could highly help conceive of strong 
interaction dynamics. 

Recognizing mass spectrum via solving non-relativistic Schrodinger equation with a 
variety of potentials is a topical issue in different branches of physics [1, 2]. It is worth 
mentioning that we focus on heavy quark mesons in solving Schrodinger equation in or¬ 
der to avoid relativistic problem. Potential models describe the mass spectrum of heavy 
quarkonium. In this system, the interaction potential is often a conhning one. 

Heavy quarkonium spectrum is a suitable instrument for exploring the static interac¬ 
tion in the infrared region. It also provides an excellent context to study the effects of 
hne and hyperhne structures. Additionally, solution of Schrodinger equation is a major 
step toward constructing a lot of potential models of quarkonium. In other words, the 
energy spectrum of the Schrodinger through various potentials is obtained in different 
ways; for example via numerical and variational methods. Hence, many potential models 
of quarkonium (such as Martin and Coulomb-Linear potential) are used that £t the mass 
spectrum of T and T [3, 4, 5, 6, 7]. 

In Sec.2 we explain our method and solve Schrodinger equation for heavy-quark mesons 
bound state. In Sec.3 introducing Be system, we obtain energy levels. Sec.4 is devoted 
to the estimating of the energy levels of ti, and in Sec.5 we present stability intervals for 
4/ and T for Yukawa plus linear potential. 


2 Lippmann-Schwinger Equation 

Since the binding energies in heavy-quark mesons are small compared to their constituent 
quark masses, one should use relativistic equations to analyze these systems. As we 
know, relativistic equations are too complicated. Hence, heavy-quark mesons (which 
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have high quark masses compared to their energy levels) are preferred to be used in the 
present study. Furthermore, Schrodinger equation, the integral form of which is Lippman- 
Schwinger equation, is employed for the calculations. 

We can use homogeneous Lippman-Schwinger equation for two particles. Schrodinger 
equation will be: 

I >= I iJb > (1) 

Equation (1) in conhguration space turns out as: 


= —m^Ti/2J dr'r'‘^J dx'J dfj)' 


f 2 n exp(—\/m\Ei,\\r — r'\) 

y I fell ^2) 


^ _ J.' 


where Ei, denotes the binding energy of the system. Equation (2) could be rewritten as: 

POO , rl 


'^bir) = dr dx M{r,r ,x )ipb{r ) 
Jo J-i 


where: 


, exp{{-Jm\Eb\)Vr^ + - 2rri) 

M{r,r ,x ) = —mv27i2n - . - V{r ) 


— 2rfx 


(3) 


(4) 


The eigen-value form of Equation (3) is: 


KiEb)]^^^ >= KEb)\'iPb> 


(5) 


where A = 1 is the positive eigen-value of the highest value to solve this equation. We 
use Gauss-Legendre method and integral equation is transformed into [8]: 


K K 

Mr) = -mV^2,rY.'£KX% 

j=li=l 


V2 


exp{-^Jm\Eb\p{r, r-,3;))) 


p(r, r'i.x'j] 


V{ri)Mrd (6) 


Equation (6) could be rewritten as: 


Mr) =iN{r,r[)iJb{r'i) (7) 

2=1 


that: 


N(rX) = 

i=i 


exp{-^Jm\Eb\p{r, r', x' )) 
p(r, r'i,x'j) 


VM 


( 8 ) 


Then we use a Fortran code that renders the kernel in Lippman-Schwinger equation di¬ 
agonal and presents the eigen-valnes as its outpnt. Appearance of A = 1 in the outpnt 
indicates that the energy supposed for the system is the binding energy. This way, we 
obtained the energy levels for tt and be. 
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3 Be Energy Levels 


Investigation of the spectrum of system is a big step toward better calculation of 
the quantitative properties of quark models and improving QCD sum rules. The B^. 
quarkonium states are capable of submitting an excellent representation of heavy quark 
dynamics. 

The B(. meson is fascinating in that it enables us to employ phenomenological information 
derived from the investigation of Bottomonium and Charmonium. 

This meson is the only heavy meson which is made up of two heavy quarks with different 
flavors of b and c. Charm and Beauty are regarded as heavy mesons because their masses 
highly exceed the scale of interaction of QCD {mg 3> ^qcd) [9]- The mass spectrum 
of this meson lies somewhere between those of ip{cc) and T{bb) systems. Therefore Be 
meson is an interesting topic for study. Be system has been studied experimentally [10] 
and elaborated on theoretically [11, 12, 13, 14, 15, 16]. 

In order to calculate the energy levels of the Be system, we use the Martin potential with 
these data [6]: 

V{r) = A{rAy + B 
rrie = 1.65 GeV 
rrib = 5.04 GeV 

A = 8.068 GeV (9) 

A = 1 GeV 
B = 6.869 GeV 
u = 0.1 


We seek A = 1 through the eigen-value spectrum in the output. The input energy in¬ 
troduced into the program that creates A = 1 is the binding energy of each level. Our 
hndings in table(l) show the energy levels of the Be system. In this table, the relativistic 
corrections are ignored. The results derived from our procedure are consistent with [17] 
and [18]. 

Table (1) shows that the results with nS states are more consistent with [17, 18] than the 
results with nP and nD states. 
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Table 1: Be Energy Levels 


state 

1^ 

2^ 

3^ 

ref[17] 

6.301 

6.893 

7.237 

ref [18] 

6.344 

6.910 

7.024 

this method 

6.28 

6.93 

7.12 

state 

2P 

3P 

AP 

ref[17] 

6.728 

7.122 

7.395 

ref[18] 

6.763 

7.160 

— 

this method 

6.74 

7.2 

7.29 

state 

3D 

AD 

5D 

ref[17] 

7.008 

7.308 

7.532 

ref[18] 

7.030 

7.365 

— 

this method 

7.19 

7.25 

7.34 


4 Predicting Topomonium Energy Levels 

A number of factors involving hyperfine splitting, spin orbit, and revising can result in 
better modeling of ip{cc) and T{bb). Of course, revising different models is not necessary 
for topomonium since we want to obtain the approximate position of energy in topomo¬ 
nium spectrum [19]. 

According to phenomenological studies, t-quark has a short half-life [20]. We consider its 
mass spectrum in this short half-life. 

In the present work we employ Coulomb-Linear potential with this form [4, 5]: 

V{r) = -^ + Br + C 
C = -0.542 GeV 

B = 0.185 GeV^ (10) 

A = 0.47 

rrit = 21 GeV[lA\ 

Cornel potential consists of two terms. One is Coulomb potential, which applies to short 
distances. This term displays the exchange of one gluon. The second term applying 
to large distances is linear potential that is indicative of quark conhnement. We solve 
Lippman-Schwinger equation to detect the eigen-value of A = 1. We estimate the energy 
levels of topomonium in table (2). 
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Table 2: mass spectrum of topomonium with Coulomb-Linear potential coefficient mt = 
2lGeV, A = 0.47, B = 0.185GeV^, C = -0M2GeV _ 


state 

1^ 

2S 

3S 

AS 

5^ 

6A 

7S 

ref [19] 

40.352 

41.372 

41.701 

41.924 

42.108 

42.269 

42.407 

this method 

40.1 

40.45 

41.9 

41.98 

42.2 

42.39 

42.8 

state 


2P 

3P 

AP 

5P 

6P 

7P 

ref [19] 


41.33 

41.699 

41.894 

42.075 

42.233 

42.376 

this method 


41.02 

41.39 

42.01 

42.19 

42.73 

42.85 

state 



3D 

AD 

5P 

6D 

7D 

ref [19] 



41.615 

41.846 

42.029 

42.188 

42.335 

this method 



41.79 

41.8 

41.10 

41.5 

42.85 


Our results are consistent the previous hndings in the literature [19]. It is worth 
mentioning that n-th state masses are obtained from the binding energy according to: 


Mn,l = 2mq + En,l 


( 11 ) 


As we know spin dependent interaction in heavy systems is small (of the m~‘^ order) 
and can be ignored. The binding energy can also be used to calculate the spin-coupling 
constant A which is dehned as: [8, 21] : 


A 


^ mim2 

S1.S2 


( 12 ) 


5 Stability Interval 

Yukawa plus linear potential is indicative of the appropriateness of our adopted method. 
Yukawa plus linear potential [1] is of this form : 

V [r) = a^r + Vi^—+ V 2 ^ 

Oi = IGeV'^, ui = 10, fii = SGeV, = IGeV, V 2 = —5 

The eigen-values of A are calculated through solving Schwinger equation. We study the 
range of r' around r-cutoff [8]. (R-cntoff results from introducing an artihcial barrier for 
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Table 3: 

stability interval of T and 4/ 

Spin 

Meson 

Eb 

Stability interval (fm) 

S=1 

cc 

-206.63 

[1.2, 2], [2.5, 5.2] 

S=1 

hi 

-891.93 

[1,2], [2, 2.5] 

s=o 

cc 

-619.91 

[1,2.1], [2.2, 5.2] 

s=o 

hi 

-573.3 

[0.22, 2], [2.2,12] 


the potential). The interval in which A = 1 will be stability interval of this potential. This 
interval shows the range where all potential-coefficients, masses, and binding energies are 
identihed accurately. We show the intervals for Charmonium and Bottomonium in table 

(3). 

6 Conclusion 

Through the procedure discussed, we solve the Schwinger equation using Martin and 
Coulomb-Linear potential for the purpose of solving identify the eigen-value problem of 
quark and antiquark bound states. With this very method, we managed to the energy 
levels of be and ti heavy mesons. The results comply with those in the literature. 
Binding energy is important in view of the fact that it could be employed to obtain the 
n-th state mass and coupling coefficient. Through the analysis of the eigen-value spec¬ 
trum, we identihed the stability intervals in which all the potential coefficients of Yukawa 
plus linear potential, masses, and binding energies are well identihed. 
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